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Abstract 

We investigate the nonlocal behavior of passive tracer dispersion with ran- 
dom stopping at various sites in fluids. This kind of dispersion processes is 
modeled by an integral partial differential equation, i.e., an advection-diffusion 
equation with a memory term. We have shown the exponential decay of the 
passive tracer concentration, under suitable conditions for the velocity field 
and the probability distribution of random stopping time. 
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1 Introduction 



The discharge of pollutants into coastal seas or rivers is common but harmful. The 
pollutants may be held up or trapped at some sites during the process of their 
dispersion. For the benefit of better environment, it is important to understand the 
dynamics of such passive tracers. 
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The Eulerian approach for studying passive tracer dispersion attempts to under- 
stand the evolution of tracer concentration profile as a continuous field quantity 



12], which ultimately, at large times, satisfies the advection-diffusion equation. 

When there are side branches, traps or some special sites in a shear flow, where 
passive tracers may be held-up or arrested and the stopping times of tracers at these 
traps are random, the tracer concentration profile C(x,t) then satisfies a nonlocal 



transport equation derived by Young |T6| 



C t + U{x)C x = DC XX - d t f k{t - s)C{x, s)ds (1) 

Jo 

where D > is the diffusion constant and U(x) is the smooth fluid velocity field. 
The kernel k(t) > is proportional to the probability distribution of tracer stopping 
times at traps. 



Young [16j has shown that the anomalous diffusion phenomenon can occur under 
certain conditions on velocity field U(x) and kernel k(t). 

In this paper, we study the long-time behavior of the concentration profile. We 
have shown the exponential decay of the passive tracer concentration, under suitable 
conditions for the velocity field and the probability distribution of random stopping 
time. 

2 Exponential Decay with 

Increasing Random Stopping Distribution 

We rewrite equation ([!]) as 

C t + U{x)C x = DC XX - k(0)C - f k'(t - s)C(x, s)ds, (2) 

Jo 

where x varies on the real line R = (—00,00). We impose an appropriate initial 
condition 

C(x,0) = C Q (x). (3) 

We denote L 2 (R) as the standard space of square- integrable functions. The usual 
mean-square norm in this space is denoted as || ■ ]]• Moreover, denotes the 

usual Sobolev space. Note that from [|TJ, p. 57, H 1 ^) = Hq(R). We also denote 
L 1 (R + ) as the space of integrable functions on positive real line R + = (0,oo). 
Moreover, C 1 (R + ) is the space of continuously differentiable (smooth) functions on 
R+. 

In order to study the evolution of the concentration C(x,t), we should estimate 



the nonlocal convolution integral in (Eh, and thus we need the following result [14 . 
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Lemma. If k 6 I}(R + ) zs a positive kernel and satisfies k' G L 1 (R + ), then for any 
yeL\R+), 

I \k * y(r)\ 2 dr < [3qK / k * y y(r)dr, 

Jt JtQ 

where K = \k\ 2 + i\k'\ 2 , < t < h < oo, and (3 > is such that k(t) — /^e - * > 0. 
Moreover, if k(t) — /^e -7 * > for some constant 7 > 0, then the conclusion also 
holds with (3qK replaced by ^PoK. 

Now we estimate the norm of the concentration C(x, t). 

Taking the inner product of ([!]) with C in L 2 (R) and integrating by patrs, we 
obtain 



1 d 1 r°° 

-j t \\Cf + D\\C X \\ 2 - - J ^ U\x)C 2 dx + k(0)\\C\\ 2 = 

/oo rt 
/ k'(t- s)C(x,s)dsC(x,t)dx. (4) 
-00 Jo 



We assume that 

-U'(x) < —a + k(0), for some positive constant a , (5) 

k' > 0,and k, k', k" e L 1 (R + ), (6) 

which imply that the velocity field U(x) has bounded gradient and the probability 
distribution of the random stopping times is strictly increasing. 

Integrating (f|) with respect to t from t to £ l5 and using @, @ and the Lemma 
above, we obtain that 

1 (iicc*oir - nc(*b)ir) + «b Z* 1 nee*)!! 2 *, < 0. (?) 



2 Jt 

Thus we have 

||C7(*i)H a < ||C7(*o)|| a , V * <ti. 
If we take t = t, ti = t + 1 in (0), then 

\\C(t + 1)|| 2 < T 1 \\C(s)\\ 2 ds < ^-(\\C(t)\\ 2 - \\C(t + l)f), (8) 



and 



max \\C(s)\\ 2 = \\C(t)\\ 2 = \\C(t)\\ 2 - \\C(t + 1)|| 2 + \\C(t + l)f 

<<s<t+l 

<(i + ^)(ll^)ll 2 -lic(t + i)|| 2 ). 



By the difference inequality of Nakao JTO] |TI[ , we finally obtain 



\\C{t)f < Ae~ B \ t>0, (9) 

for some positive constants A and B. Thus ||C(£)|| decays exponentially fast. There- 
fore we have the following result. 

Theorem 1 If the velocity field U(x) has bounded gradient and the probability dis- 
tribution k(t) (of the passive tracer random stopping times) is strictly increasing, 
namely, they satisfy conditions ^) and respectively, then the concentration of 
the passive tracer with random stopping approaches zero exponentially fast in mean- 
square norm. 

3 Exponential Decay with History Source Term 

If there is a source term due to the past history of C from — oo to 0, more precisely, 
if there is a source term 

r° 

f(x, t) = - / k(t - s)C(x, s)dx, x E R, t > 0, 

J — oo 

then, d||) can be written as 

C t + U{x)C x = DC XX - k(0)C - f k'(t - s)C(x, s)ds 

POO 

= DC xx -k(0)C- k'(s)C(x,t- s)ds. (10) 
Jo 



Let 



C t (x,s) = C(x,t-s), (11) 
n t {x,s)= f 'c t {x,T)dr= [ C{x,r)dr, s>0. (12) 

JO Jt-s 

We further assume that lim k'(s) = 0, then 



s— >oo 



k'(s)C(x,t — s)ds = / k'(s)C t (x, s)ds = — k"(s)i] t (s)ds. 

o JO JO 

Let /i(s) = k"(s), then (|1^) can be reformulated as 

roo 

C t + U(x)C x = DC XX - k(0)C - / /i(a)j7*(ar, s)ds, (13) 



o 



rf t (x,s) = C{x,t)-^ s rf. (14) 
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The initial conditions are then 



C(x, 0) = Cq(x), rj (x, s) = rjo(x, s) = / C(x,t)cIt, x E R. (15) 

J — s 

Let 

H = L 2 {R) x L 2 {R + - //; L 2 (R)), 

where 

L 2 (R + ; fx- L 2 (R)) = {ip:R + ^ L 2 (R)\ / n{s)\\ip{s)\\ 2 ds < oo}, 

Jo 

and < ■ > fi , || • \\fj, denote the inner product and norm in this space, respectively. 
The well-posedness of the initial boundary value problem of ([13]), (|14]) and (|15|) 



can be shown as in [0] and ||. Here, we focus on the asymptotic behavior under 
the condition (|5|) on the velocity field U(x), and the following conditions on the 
probability distribution k(t) of the random stopping time 

k' > and lim k'(s) = 0, (16) 

k" E C 1 ^) flL 1 ^) and k"{s) < 0,k"'(s) > 0, Vs G (17) 
fc"'(s) + <yjfc"(s) < 0, Vs G and for some constant 5 > 0. (18) 

We now estimate the evolution of C, 77* described by (|T3|)-(|l4|). Multiplying (p~3|) 
by C and ([14]) by 77', then integrating by parts and adding the results, we get 

~(ll^l| 8 + 11^11;)-/^^)!^ 

+ fc(0)||C|| 2 + D\\C X \\ 2 + < ^-rf,r? >,= 0. (19) 

Since for a fixed t G (0, 00) both /u||t7*|| and Jj»7*|| belong to L 1 (R + ), we now have 

hmMWOOir = 0. (20) 
Hence, using ( P0|) and integrating by parts, we obtain 

<|^W>m=-/ vt(*W\\ 2 ds. (21) 
By putting (§), QI§)-(]ig), and (0]) into (0), we finally get 

|(I|C|| 2 + \\V%) + 2«o||C|| 2 + 25||t/|| 2 < 0, (22) 



where a , 5 are defined in (0), fll8]) , respectively So by the Gronwall inequality 
we obtain the exponential decay for ||C|| 2 + 1 1 77* || 2 and thus for the concentration (in 
mean-square norm) ||C||. We summarize this result in the following theorem. 
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Theorem 2 Assume that the velocity field U(x) has bounded gradient and the prob- 
ability distribution k(t) (of the passive tracer random stopping times) is strictly 
increasing, namely, they satisfy conditions ffi) and ^), respectively. If k(t) satisfies 
further growth conditions flTb)), ^Tlj ) and (fiTj), then the concentration of the passive 
tracer with random stopping and with a history source approaches zero exponentially 
fast in mean-square norm. 

Note that the transformation that leads to the reformulation (0)-(0) is inspired 
by H . Under this transformation, a nonautonomous system ( ]T0| ) with history source 
term becomes an autonomous sysytem (|13|)-(|14l), at the expense of increasing system 
components or dimension. 

4 Discussions 

In this paper, we have studied the nonlocal behavior of passive tracer dispersion with 
random stopping, using an advection-diffusion equation with a memory or history 
term. We have shown the exponential decay of the passive tracer concentration, 
under suitable conditions for the velocity field and the probability distribution of 
random stopping time. There are other situations where tracer dispersion is modeled 
by advection-diffusion type of equations with memory terms JHJ, || and Q. 
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